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With an effective Lagrangian approach, we give a full analysis on the NN — > NNtttt 
and NN — > NNtttt reactions by exploring the roles of various resonances with mass 
up to 1.72 GeV. We find large contributions from A, AT* (1440), A(1600) and A(1620) 
resonances. Our calculations also indicate sizeable contributions from nucleon poles for 
the energies close to the threshold. A good description to the existing data of different 
isospin channels of NN — ¥ NNtttt and NN — > NNtttt for beam energies up to 2.2 GeV is 
reached. Our results provide important implications to the ABC effect and guildlines to 
the future experimental projects at COSY, HADES and HIRFL-CSR. We point out that 
the PANDA at FAIR could be an essential place for studying the properties of baryon 
resonances and the data with baryon and anti-baryon in final states are worth analyzing. 
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1. Brief summary of current status 

Our understanding on the low energy strong interaction physics is still unsatisfac- 
tory though quantum chromodynamics (QCD) has been established as the standard 
theory of strong interaction for many years. Because of its color confinement non- 
perturbative properties, we have difficulties to calculate the properties of mesons 
and baryons directly from QCD. Constituent quark models are developed to ex- 
plain the mesonic and baryonic spectrum and they are successful in some aspects. 
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But one of severe problems is that they predict more excited states around 2 GeV 
than what have been observed experimentally. As a result, a lot of efforts have been 
devoted to the meson production in pion-, photo-, and electro-induced reactions in 
order to reduce the uncertainties of extracted parameters of resonances, deepen our 
knowledge on the structure of the resonances and search for missing resonances. 

In the past few years, the double pion production in nucleon-nucleon collisions 
tends to be a fascinating field for studying resonances properties and has been 
accurately measured at the facilities of CELSIUS^, COS^, KEEP, and PNPI- 
Gatchina^. In the TableJTJ we list all the data measured after the year of 2000. The 
CELSIUS^ and COS v2l Collaboration measured the differential cross sections of 
pp — > ppir + ir~ and pp — > ppir°ir° channels for the beam energies from the thresh- 
old to 1.3 GeV. One of the promising findings was that the iV*(1440) contribution 
dominates at the close-to-threshold region, as expected by the Valencia modeM 
This provided us a good place to explore the properties of iV*(1440) whose theoret- 
ical interpretation is still under controversial. However the Valencia modeL though 
compatible with the old bubble chamber and magnetic spectrometer data^, over- 
estimated the new data at the close-to-threshold energies by several times. Another 
finding was that the CELSIUS data of irir invariant mass spectra for the beam ener- 
gies above 1.0 GeV demonstrated a single peak at the low invariant mass while the 
model gave double hump structure which was inconsistent to the data. The model 
also predicted a preferential parallel emission of the two pion-meson, which was 
contrary to the CELSIUS data. The third interesting finding was that there was a 
level-off behavior around the beam energies of 1.0 GeV in the energy dependence 
of the pp — > ppir°ir°, which was thought to be the result of interference between 
different contributions. KEK Collaboration made another progress and they mea- 
sured the total cross section of pn — > prnx^ix^ and pn — > ppir^n channels which 
were important for our understanding on the excitation of AT* (1440). 



Table 1. The data of nucleon-nucleon collisions measured after the year of 2000. Those with 
the data of differential cross sections are marked by bold characters for the beam energies. 





Channel 


Collaboration (Tp (MeV)) 


PP - 


-> ppTT + 1T~ 


CELSIUS(650, 680, 750, 775, 895, 1100, 1360), 






Gatchina(717, 818, 861, 900, 980), COSY(750, 800), 






KEK(698, 780, 814, 908, 995, 1083, 1172) 


pp 


-> pp7r°7r° 


CELSIUS(650, 725, 750, 775, 895, 1000, 1100, 1200, 1300, 1360) 


pp - 


■+ nmz + TT + 


CELSIUS(800, 1100) 


pp —> pn7T+7r° 


CELSIUS(725, 750, 775, 1100) 


pn - 


-> pnir+n 


KEK(698, 780, 814, 908, 995, 1083, 1172) 


pn 


-> ppn~n° 


KEK(698, 780, 814, 908, 995, 1083, 1172) 



On the other side, though the role of iV*(1440) in nucleon-nucleon collisions 
was firmly established, its contribution in NN — > NNttt: reactions has never been 
considered. JETSET^l measured the pp ppir + ir channel in order to search for 



September 7, 2010 0:8 WSPC/INSTRUCTION FILE 



cx-ijmpa 



Double pion production in NN and NN collisions 3 



narrow resonances decaying to pp but they only included the double-A diagram 
in their Monte-Carlo simulation. Moreover, it has never been explored whether it 
was possible to extract the properties of other resonances from antinucleon-nucleon 
collisions. 

The one-pion exchange (OPE) models-^ of more than 45 years ago only con- 
sidered the double-A diagram and the Valencia model of more than 10 years ago 
made a further step by including also the TV* (1440) resonance and the a- and p- 
meson exchange. With the new accurate data, it is very necessary to perform a 
more comprehensive analysis including more resonances and matching all the data 
of (anti)nucleon-nucleon collisions. 

2. Description of our model 

Our full model is demonstrated in Ref. HJ so herein we only give a brief description 
to the main features. The effective Lagrangians for the reson ances in our model are 
based on a Lorentz covariant orbital-spin (L-S) scheme In view of the overall 
system invariant mass about 2.8 GeV for Tp =2.2 GeV, we have checked contri- 
butions from all the well-established N* and A* resonances below 1.72 GeV. The 
coupling constants appearing in relevant resonances are determined by the empir- 
ical partial decay width of the resonances taken from Particle Data Group (PDG) 



A* (1620) resonances play the relatively significant role in the considered energies 
and other resonances give negligible contributions^. 

In our model, the adjustable parameters are the cut-off values in form factors at 
vertices and resonances. We include five types of form factors: those at the meson- 
(anti)nucleon-(anti)nucleon vertices, those at the meson- AT* (A*)-(anti)nuclcon ver- 
tices, the Blatt-Weisskopf barrier factors at the iV*(1440)-A-7r vertex the form 
factor at the a-n-n vertex, and the form factor of resonances and nucleon poles. 
The pp — > nnir + Tr + channel is very useful to pin down these cut-off values in form 
factors of the relevant A and A* resonances because it has negligible N* contri- 
bution, and then we could easily determine the cut-off values of N* resonances in 
other channels. 

3. Results and discussion 

If the doublc-A mechanism is dominant at high energies as in the OPE model and 
Valencia model, the total cross section of pp —> pptt°tt should be a factor of about 
four larger than that of pp — > nnn + TT + according to isospin coefficients. However, 
the new exclusive measurements indicate an approximately eq ual value of these two 
channels, which is consistent to the old bubble-chamber dataP^. After calculating 
all the possible contributions, our model got two important conclusions: 




A(1232), A*(1600), and 



(i) Our model reduces the relative branching ratio of iV*(1440) — > An and assumes 
a smaller cutoff parameter for the wNA coupling so the relative contribution 
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M„v(GeV) cos5„v 



Fig. 1. The invariant mass spectra and angular distributions of the opening angle for 7r°7r° 
systerJSl in the overall center-of-mass system for pp — > pp7r°7r° at the beam energies of 795 MeV 
(a, b) and 1300 MeV (c, d). The dashed and solid curves correspond to the phase space and full 
model distributions , respectively. In (a, b), the dotted curves correspond to TV* (1440) — > Na 
contributions. In (c, d), the dotted curves correspond to double- A contributions. 




Tp (GeV) 



Fig. 2. Total cross sections of pn — > pnw n reaction. The dotted, dash-dotted, dashed, dash-dot- 
dotted, solid and bold solid curves correspond to contribution from double- A, TV* (1440) — > Na, 
TV* (1440) — >• At, nucleon poles, full model results without and with FSI, respectively. 
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from the TV* (1440) — > Act term increases significantly, 
(ii) our model introduces significant contributions from A — > Nn — > Nnn at en- 
ergies near threshold and from A* (1600) and A* (1620) at energies above 1.5 
GeV. This is more obvious in the pp — > nnn + n + channel because the isospin 
coefficients of these terms are much bigger than those in other channels. 

As a result, our description of the total cross sections of all isospin channels is 
considerably improved and our results agree well with the new data of NN collisions 
in the close-to-threshold region. However, we find that it is difficult to interpret the 
nn invariant mass spectra for the beam energies above 1.0 GeV and the level-off 
behavior around 1.0 GeV in the total cross section of the pp —> ppTr Q ir° channel. In 
Fig. [TJ we give the invariant mass spectra and angular distribution of the opening 
angle for 7r°7r° system in the overall center-of-mass system in the pp — ¥ ppn°n° 
channel at the beam energies of 795 MeV and 1300 MeV. At 795 MeV, which is close 
to the threshold, the A* (1440) is dominant and the given spectra is consistent to the 
data. At 1300 MeV the model fails to describe the data because the A* (1440) -> An, 
which is significant at high energies, gives a double hump structure in the 7r°7r° 
invariant mass spectrum. But if the A* (1440) contribution is negligible which is 
caused by some kind of destructively interference in this energy region, then the 
double-A diagram, which is dominant by 7r-meson exchange in our model, describe 
the data welW The situation is similar in the pp — > ppn + n~ channel. 

Our model provide essential hints to the ABC effect in pn — > dn°n° reactiorfSl 
In Fig [5] we show the total cross section of pn — > pnn°n° channel. It can be seen 
that the A* (1440) -> Act and A* (1440) -> An are significant in all the considered 
energies. At the close-to-threshold region, the nucleon poles are very important. So 
it is suggested that the A* (1440) resonance and nucleon poles should be carefully 
considered before we reach the right conclusiorllSl. 

The agreement between our model and the data of NN collisions is goocPt It 
should be addressed that the A* (1440) and other resonances should be included 
in the Monte-Carlo simulation of the experimental analysis in order to get the 
correct total cross sections from limited measured phase space, which is overlooked 
by the previous measurements. We find that the pn — > pfin~n~ reaction is useful to 
determine the model parameters because like the pp — > nnn + n + reaction, the A* 
contribution is small. We also point out that other channels in antinucleon-nucleon 
collisions could serve as a good place for studying baryon resonances. 

A lot of precise measurements in nucleon-nucleon collisions are being carried out 
by CELSIUS and COSY collaborations in order to study A* and A* resonances. 
Besides, the HADES collaboration has measured some channels of AA collisions at 
the beam energies of 1.25 Ge#3 and their results are expected to come soon, so 
this will constitute a good test of our model. Recently, a cooler storage ring HIRFL- 
CSR, which can produce the proton beam with similar beam energy range of COSY, 
has already been successfully installed at Lanzhou. With the scheduled An hadronic 
detector for complete measurement of diverse differential cross sections, it will have 
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a special advantage for studying excited nucleon states through nucleon-nucleon 
collisions. We also suggest that the double pion production in antinucleon-nucleon 
collisions should be measured at PANDA (anti-Proton ANnihilation at DArmstadt) 
at the GSI Facility of Antiproton and Ion Research (FAIR) , which could provide the 
antiproton beam of kinetic energy ranging from 1 to 15 GeV with the luminosity 
of about 10 cm s . We address that PANDA could play an important role in 
the hadronic physics with a An solid angle detector with good particle identification 
for charged particles and photons. The planned new experiments are anticipated to 
offer more tools and information to help us understanding the low energy physics 
better. 
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